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A NOTE ON THE CRACKED PLATES REINFORCED BY 
A LINE STIFFENER^ 

O.S. Yahsi and F. Erdogan 
Lehigh University, Bethlehem, PA 18015 

1 . Introduction 

In this note the problem of a cracked plate reinforced by a line 
stiffener is reconsidered. The original solution of the problem was given 
in [1]. A variation of the problem with debonding between the plate and 
the stiffener near the cracked region was considered in [2]. However, 
the special case of the problem in which the crack tip terminates at the 
stiffener does not appear to have been studied.- In practice the solution 
may be necessary in order to assess the crack arrest effectiveness of the 
stiffener. In this note the problem described in Figure 1 is reformulated, 
the asymptotic stress state near the crack tip terminating at the stiffener 
is examined, and numerical results are given for various stiffness constants. 

2. The Integral Equation 

Consider the plane elasticity problem shown in Figure 1. It is 
assumed that the elastic plane is reinforced by a stiffener which is em- 
bedded into the medium. Let the thickness of the stiffener in y-di recti on 
be sufficiently small so that its in-plane bending stiffness may be 
neglected. Hence, the stiffener may be approximated by a membrane. The 
Airy stress function for the stiffened elastic plane may be expressed as 
follows: 
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where R-j , R£ and S^,...,S 4 are unknown functions. From Figure 1 one may 
note that if the loading is symmetric, x=0 is a plane of symmetry (which 
is assumed in (1)) and the six unknowns of the problem may be determined 
from the following continuity, equilibrium and boundary conditions: 


u(x,+0) - 

u(x,-0) = 

0 , 

0 <_x <°° , 

(2) 
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» yy (x.-H3) 

- o (x.-0) 
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(4) 

E C a xy( X » +0 ) - a X y( X ’- 

0)] 
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< x< °o , (5) 

a xy (0,y) 

= 0 , < 

; y < co 
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(6) 

a xx (0 >*> 

= f(y) , 

-d <y 

<-b , ] 


u(0,y) = 

0 , -°° < y 

< -d 

, -b <y <=° . j 

(7a, b) 


Equation (5) is obtained from the equilibrium of the stiffener (of length 
dx) in x direction by observing that the strain in the stiffener is 3u/3x. 
The stiffness constant y is, therefore, given by 

Y = E S A $ /E (8) 
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where E $ is the Young's modulus and A $ is the cross-sectional area of 
the stiffener corresponding to the unit plate thickness and E is the Young's 
modulus of the plane. Using the expressions 

a xx = » a yy = 3 2 4>/3x 2 , = - 3 2 <J>/3x3y (9) 

and the Hooke's Law, five of the unknowns in (1) may be eliminated by 

the homogeneous conditions (2)-(6) and the sixth may be determined from 

the mixed boundary conditions (7). 

Defining a new unknown function 

— u(x,+0) = g(y) , -°° < y < 00 (10) 

and from (7b) observing that g(y) = 0 for -®<y<-d and -b<y<«>, after 
some simple manipulations, (7a) may be reduced to 

1 f b sll^t + 1 f b k ( y}t ) g ( t ) dt = | f(y)j _ d<y< _ b f (H ) 

77 J-d 1 y 77 J-d L 

where 

r ftP ( t+ y)s i+M 

k(y,t) = 2y j Tg (it,)(3-v)+4 (1 +L T tg ) [3+v»(l»v>)yB3dB • (12) 

If the crack is embedded in the left-hand plane (i.e., if -d < -b < 0) , then 
the kernel k(y,t) is bounded for all values of y and t in the closed 
interval [-d,-b]. However, if the crack tip terminates at the stiffener 
(i.e., if b=0, or c=a. Figure 1), then it may be shown that k(y,t) becomes 
unbounded as y and t go to the end point b=0 simultaneously. The singular 
part of k(y,t) may easily be separated by examining the behavior of the 
integrand in (12) for B-x=°. Thus, it can be shown that 

k(y»t) = k s (y,t) + k f (y,t) , (13) 
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k $ (y,t) 


(3+vHl-v) _J 3(1 +v) ,y , 2(1 +v) ,y 2 

(1+v) (3-v) t+y 3-v (t+y) 2 3-v (t + y) 2 ’ 


k f (y,t) 


8 

3-v 



(t+y) 8 


[l+tg(l +v)/2] 

y(1+v)(3-v)B+4 


dB . 


(15) 


In the limiting case of b=0, k s (y,t) and the Cauchy kernel in (11) con- 
stitute a generalized Cauchy kernel. 

From (7b) and the definition (10) it is seen that the integral 
equation (11) must be solved under the following single-valuedness condition: 

f-b 

g(t)dt = 0 . (is) 

J -d 


It should be noted that equations (11 )-(15) are obtained under 

generalized plane stress assumption. For plane strain E and v should be 

2 

replaced by E/(l-v ) and v/(l-v), respectively. 

For -b < 0 the solution of (11) may be obtained by introducing the 
following dimensionless quantities: 


n 



d+b\ 

2 J • 



d+b\ 
2 ' 


9 


G(t) * g(t) , (-1 < (n,x) < 1 ) , 


(17) 


and by using a standard Gauss-Chebyshev integration procedure [3]. In 
this case the unknown function G(t) is of the following form 

G(t) = h(t)/ /- 1 2 , (18) 

where h(r) is a bounded function. After determining h(r), the Mode I 
stress intensity factors at the crack tips may be defined by and calcu- 
lated from 
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k,(-b) = Urn v£(y+b) cr (0,y) , (y>-b) , 

1 y->-b xx 

= -| Tim ^(-b-y) -J- u(+0,y) , (y<-b) , 

2 y—b 9y 

= -§ /('d-b)72 h(l) ; (19) 

k,(-d) = lim /2(-y-d ) a (0,y) , (y<-d) , 

1 y-*-- d xx 

= | lim /2(y+d) u(+0,y) , (y>-d) , 

y-*-d ojr 

= | /(d-b)/2 h(-l) . (20) 

In the case of b=0, because of the kernel k $ given by (14) the solution 
of the integral equation (11) has no longer square root singularities. 

By introducing the new normalized quantities 

n = f y + l » T “ § t + 1 , G(t) = g(t) , -1 < (n.x) <1 , (21) 

assuming G(i) of the form 

G(x) = h(x)(l-T) a (l+r)^ , (-1 < Re(a,S) <0) , (22) 

and by following the function theoretic method outlined, for example, 
in [3] the characteristic equations of the problem giving a and p may be 
obtained as follows: 

cotTrS = 0 , (23) 

cosTia + ^ (a 2 +2a) - [^j{{ = 0 . (24) 

Equation (23) gives the expected result of 3 =-1/2 whereas (24) shows that 
a depends on the Poisson's ratio of the plane only. Again note that (24) 
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is obtained for plane stress; for the plane strain case v should be 
replaced by v/(l-v). With a and B known, (11) may be solved by using a 
Gauss-Jacobi integration formula [3], 


3- Asymptotic Stress Field 

For the plane with a crack terminating at the stiffener in order to 
study the initiation of various modes of fracture growth the asymptotic 
behavior of the stress state around the crack tip x=0, y=0 has to be 
investigated. Going back to the basic formulation of the problem various 
stress components around the crack tip may easily be expressed in terms 
of bounded integrals with g(y) defined by (10) as the density function. 
For example, after some relatively simple manipulations the stress 


a xx (0,y) for y>0 may be expressed as 


f (° r 

* Jf/Jt 


4v _L 3(l+v) 

j){ 3-v) t-y 3-v 




:7T( d-V J 


:ir( 3-v 


[ d h 1+ ^-(l-) 7^7‘^ r+(3+V) t-y _ 9(t)dt 

(my - j (l - tje -rn ^ t-y ^Ei (m(t-y))g(t)dt , 


where Ei(x) is the exponential integral and 


0 <y <« , 


If we now observe that (see (21) and (22)) 

g(t) = h(T)(l-t)°()+T)*‘ 1 = G(t)(-t) a (t+d)-' i , (0 > a > -1 ) 
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and define the sectionally holomorphic function 


F<z> = ? C, w dt ■ < 28 > 

we obtain the asymptotic behavior of F(z) near the end point z=0 as follows: 

F(z) = + F (z) (29) 

v ' ^ si nuct o v ' 


where F q (z) is bounded near z=0. From (28) it is seen that at z=y>0 
F(z) is holomorphic and 


F(y) 


1 

"IT 


f° sill 
-d 


dt 


6(0) y a 
sinua 


+ 


(30) 


From (25) and (30) it can now be shown that only the terms in the first 
integral in (25) contribute to the stress singularity and by using (30) 
these terms may be expressed as 


I 

IT 

1 

TT 



— y 2 g(t)dt 
(t-y) 2 

2 

y 3 g(t)dt 
(t-y) 3 


y jL p(y\ = a Q.(Q) )L + p (y) 

y dy sinma h * 


(31) 


1 ^ p(y) = a ( a ,~T ). ^(0) y + p (y) (22) 

2 dy 2 2 ft sinTO 2 ’ } 


where F-j and ?2 are bounded near y=0. 

By substituting from ( 30)- ( 32 ) into (25), the leading term of a ^ 
near y=0 is then found to be 


a xx (0 >*> 


G(0)£ 


2/d sirnra *- 


4v 


1+v 


.1 




(33) 


where 6(0) is given in terms of the bounded function h(x) as (see (27)) 


6(0) = (2/d) a ‘ % h(l) . 


(34) 
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Similarly, other critical stress components and their asymptotic 

values near the crack tip x=0, y=0 may be determined as follows: 

C fO r /I 1 n/1. \ . 2 x. 


a yy^ 0,y ) " 2ir(3-v] 


f J4 _L_iI+vbl. 

J-dl 1+v ^ (t-y) 2 


(t-y) (t-y) 


m(l+v)t\ 

(t-y) 2 


- *4^* - - 2 »(kH(' -) * 

* Ei(m(t-y))|g(t)dt , 0<y<», 


Oyy(O.y) « 
yy 2v^ s 


Ml 
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1+v 2a n 

3-vT ~ 3^v a y 1 y>0 ’ 




TT( J-V) 


-dlt 2 + x 2 (t 2 + x 2 ) 2 


m(l-vj 

4 


H 2 + x 2 


+ ("rt-if^) £ eBt g C -^ SBx <«]j g(t)dt , 0 <x <= , 


°yy (x ’-°> 


<? X y( X »“0) 


- EG (°>, lia x° , x>o , 

2 IS sin f 3 ‘ v 


:tt( d-v 


3 j 7+2v-v 2 x 2(5-v)xt 2 2mxt 
-dl x 2 n 2 ' (x 2 +t 2 ) 2 ' tV 


+ [2m 2 t -4m(l+v)] — - g S ^ nBx dgj g(t)dt , 0<x<=°, (39) 


r (x,-0) s M^r_J5 I v]a_ + 7+2^ 

Xy 4/J L(3-v)sin ^ (3-v)(l+v)cos ^ 


! x , x>0 . 


From the asymptotic expressions (33), (36), (38) and (40) it is seen 
that stress components near the crack tip has the form 


5 r* fjjfv) = Eh(l )2 a-1 ' s (r/d) a f, j(v) , (i,j = x,y), (41) 
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where r is a small distance from the crack tip, and f . . are known 

* J 

functions and -1 <Re(a) <0. Thus, the calculated quantity h(l) is the 
measure of the stress singularity near the crack tip. For b=0 in the 
numerical analysis, in addition to k^(-d) defined by (20), the following 
"stress intensity factor" is defined andtabulated (Figure 1) 


k-| (0) 


lim SI y“ a 
y-H-0 


o xx (°.y) 


am 

_ a • 

2a sinira 


4v 

jl+v)(3-v) 


, 1 +V ( f. 

+ 3^7 (a 


+ 2a) 




(42) 


We note that for any external load leading to the opening of the 

crack G(0) <0. Since a is also negative, from (33), (36), (38) and (40) 

it follows that near the crack tip we have a (0,y) >0, a (0,y) >0, 

xx yy 

a (x,-0) >0, and cr (x,-0) <0, which are the intuitively expected results. 

XX Xjr 

The power of singularity a defined by (27) and calculated from (24) 
is given in Table 1. The stress intensity factors calculated for constant 
crack surface pressure cr (0,y) = -a , (-d<y<-b), v = 0.3, generalized plane 

XA w 

stress, various crack locations, and various values of the stiffness con- 
stant y are given in Table 2. Note that, as expected, the stress intensity 
factors decrease as y increases and as the crack moves closer to the 
stiffener. 
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Table 1. Variation of the power a of the 
stress singularity with Poisson's ratio. 



a 


Plane Stress 

Plane Strain 

0 

0 

0 

0.1 

-0.1026 

-0.1096 

0.2 

-0.1547 

-0.1743 

0.3 

-0.1913 

-0.2267 

0.4 

-0.2197 

-0.2739 

0.5 

-0.2429 

-0.3196 
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Table 2. The calculated stress intensity factors: k' (-b) = (-b)/c 0 ^a 

k'(-d) =k 1 (-d)/o 0 v€, k'(0) =k 1 (0)/a o a" a , a = (d-b)/2, c-(d+b)/2, y = E s A s /E, 
v = 0.3, ct = -0.1913, the case of generalized plane stress. 



k'(0) 


k'(-b) 




k'(-d) 



il 

1 

1.1 

1.5 

2.0 

10 

1 

1.1 

1.5 

2.0 

10 

0.1 

1.197 

0.952 

0.993 

0.998 

1.000 

0.988 

0.994 

0.998 

0.999 

1.000 

0.2 

1.008 

0.919 

0.988 

0.996 

1.000 

0.981 

0.989 

0.996 

0.998 

1 .000 

0.25 

0.953 

0.905 

0.985 

0.995 

1.000 

0.977 

0.986 

0.995 

0.998 

1.000 

0.3 

0.911 

0.894 

0.983 

0.995 

1.000 

0.975 

0.984 

0.994 

0.997 

1.000 

0.4 

0.850 

0.874 

0.978 

0.993 

1.000 

0.970 

0.981 

0.993 

0.997 

1 .000 

0.5 

0.806 

0.858 

0.974 

0.991 

1.000 

0.966 

0.977 

0.992 

0.996 

1.000 

0.6 

0.774 

0.844 

0.971 

0.990 

1.000 

0.962 

0.974 

0.990 

0.995 

1.000 

0.7 

0.748 

0.833 

0.967 

0.989 

1.000 

0.959 

0.972 

0.989 

0.995 

1.000 

0.8 

0.726 

0.823 

0.964 

0.988 

1.000 

0.956 

0.969 

0.988 

0.994 

1.000 

1.0 

0.694 

0.806 

0.959 

0.986 

1.000 

0.951 

0.965 

0.985 

0.993 

1.000 

1.25 

0.664 

0.790 

0.953 

0.983 

1.000 

0.946 

0.960 

0.983 

0.992 

1.000 

1.5 

0.643 

0.777 

0.948 

0.981 

1.000 

0.941 

0.957 

0.981 

0.991 

1.000 

1.75 

0.626 

0.766 

0.944 

0.979 

1.000 

0.938 

0.953 

0.979 

0.990 

1.000 

2 

0.613 

0.757 

0.941 

0.978 

1.000 

0.935 

0.951 

0.978 

0.989 

1.000 

2.5 

0.593 

0.744 

0.933 

0.975 

1.000 

0.930 

0.946 

0.975 

0.987 

1 .000 

3 

0.578 

0.733 

0.930 

0.973 

1.000 

0.925 

0.942 

0.972 

0.986 

1.000 

3.5 

0.567 

0.725 

0.926 

0.971 

1.000 

0.923 

0.939 

0.971 

0.985 

1.000 

4 

0.558 

0.718 

0.922 

0.969 

1.000 

0.919 

0.936 

0.969 

0.984 

1.000 

5 

0.545 

0.708 

0.917 

0.966 

1.000 

0.915 

0.932 

0.966 

0.982 

1.000 

7 

0.529 

0.695 

0.909 

0.962 

1.000 

0.909 

0.926 

0.962 

0.980 

1.000 

10 

0.516 

0.683 

0.902 

0.958 

0.999 

0.904 

0.921 

0.958 

0.977 

0.999 
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Figure 1. The geometry of a stiffened plate containing a crack. 
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